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ABSTRACT 


Dynamics of the circulation in the Sea of Marmara are 
investigated with a time dependent, three dimensional 
numerical model. The empirically inferred hydrologic regimes 
of the sea and connective straits are discussed. 

A baroclinic model based on the primitive equations is 
solved by direct integration of an initial value problem. 
The circulation in the Sea is driven by surface forces that 
Simulate wind stress and horizontal pressure eradiene forces 
related to internal stratification. 

The predicted density field, in sigma—t unitS,1is compared 
With data. Detailed three dimensional horizontal velocity 
patterns and vertical velocity patterns in horizontal planes 
are given. 

Bottom friction, irregular bottom topography, nonlinear 
terms in the momentum equation and vertical mean part of the 
horizontal velocity have been omitted. For simplicity 


density is predicted in place of temperature and salinity. 
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I. INTRODUCTION 


There is an intense growing interest in recent years in 
the development of simulation of hydrological conditions in 
coastal waters and smaller adjacent seas. Advancement in 
modern technology requires much greater accuracy than can be 
obtained by the use of classical instrumentation and measure— 
ment techniques only. Due to the discontinuous nature of 
such observations these measurements can provide the basis 
for analysis of the current system, temperature, salinity 
and other variables on a climatological basis. But climato— 
logical conditions are not accurate enough in coastal waters 
for uses such as construction, navigation, tracing of pollu- 
tants, mining effluent, oil spills, search and rescue opera- 
tions, agriculture, etc. It is also true that sampling the 
ocean adequately whether for exploratory or research purpose 
is a major problem. Great short—term variability of condi- 
tions in coastal waters makes it even more difficult to make 
adequate representative measurements. 

As a result, in addition to the classical methods of 
analyses of current systems, determination of salinity and 
temperature distributions etc., new simulation techniques 
Which were applied long before synoptic numerical weather 
prediction received much more attention (T. Levastu, S. 


Larson and K. Rabe, 1976). 
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The objective of this study is to develop a regional 
circulation model of the Sea of Marmara. The circulation 
and density distribution, in sigma—t units, in that sea are 
Simulated by a nine~level numerical model with an hori- 
zontal flat bottom at 1000 m depth. The motion is driven 
by prescribed wind stress and horizontal pressure gradient 


forces due to internal stratification. 
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II. PHYSICAL GEOGRAPHY AND HYDROLOGY OF 


THE SEA OF MARMARA 


The Sea of Marmara iS a small sea located between Asia 
Minor and Europe and occupies the main part of the Turkish 
Straits system (210 km). The Strait of Canakkale (Darda- 
nelles) (60 km) connects the Sea of Marmara to the Mediter- 
ranean and through the Strait of Istanbul (Bosporus) (30 km) 
it is connected to the Black Sea at the North. (Figure 1) 

The structure of the baSin which is small in area (about 
11000 km} is associated with the Anatolicin fault which 
lies along its bottom to the North. The trough passing along 
the Northern deep slope consists of three depressions (max. 
1380 m). These three baSins are separated by low connecting 
Sills. A shallow trough lies at the foot of the continental 
Slope which is well pronounced only in the easternmost part 
of the sea. At the South it becomes sHallower having depths 
of the order 60-80 m. Minimum depths at the Strait of 
Canakkale and Strait of Istanbul are 57 and 37 m; they have 
mean widths 4.5 km and 0.7 km and average lengths 60 and 30 
km, respectively. 

The climate of the Sea of Marmara is influenced by the 
regional atmospheric circulation system of the Eastern Medi-— 
terranean, the Black Sea and Asia Minor Peninsula. During 


the summer, the Eastern Mediterranean basin is under the 
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influence of a stable wind system which is called "Etesian” 
wind. EteSians are monsoon~like and associated with a quasi- 
Stationary cyclonic circulation above the Asia Minor Penin- 
Sula during summer. The active period of the Etesian winds 
falls almost entirely within the summer season (approximate- 
ly from May to October). There are two maximums, which are 
pronounced during July and August. These Etesian winds, 
which are dry winds and generally northerly or northeasterly 
from the Black Sea, have diurnal variability. 

Southerly winds bring warm and humid air to the area 
during winter. When these southerly winds are sufficiently 
Strong, they produce storm waves and surges along the coast 
of Marmara and at the entrance to the Strait of Istanbul 
(Gunnerson and Ozturgut, 1974). 

The hydrological regime of the Sea of Marmara is largely 
determined by the water exchange between the Mediterranean 
and Black Seas through the connective straits, Strait of 
Canakkale, and Strait of Istanbul. Problems of water ex 
change, vertical mixing and hydrological characteristics of 
the currents in these straits and in the Sea of Marmara have 
been the subject of many extensive studies conducted by the 
Turkish Navy Hydrographic Office and by the University of 
Istanbul. In addition to these fundamental and frequent 
Studies there were some investigations conducted by German, 
Russian and English investigators at the end of the 19th 
and beginning of the 20th century. Recent studies and obser— 


vations, if connected to old data and obServations, can give 
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a better understanding of the complete regional oceanographic 
condition. 

These Studies are mainly concentrated on the Strait of 
Istanbul and the vicinity of this strait. This is due to the 
importance of this strait as a link in the connection between 
the Black Sea and the Mediterranean Sea through the Sea of 
Marmara. This makes it an important factor controlling Black 
Sea chemistry, biology, sedimentation, etc. 

Based on the first detailed survey of Makarov and Merz 
1917-18, in Moller (1928), and subsequent studies of Ullyott 
and Ilgaz (1946), Pektas (1956), A. K. Bogdanova (1961, 1965), 
C. G. Gunnerson and E. Ozturgut (1974) and others, the main 
properties of this straits! system and its vicinities are 
fairly well understood. 

One of the main characteristics of these straits isa 
two—layer structure of the current, which is determined by 
the differences between the relatively fresh water of the 
Black Sea and saline water from the Mediterranean which occu- 
pies the depths of the Sea of Marmara and forms the bottom 
water of two connective straits (Figure 2). Other external 
factors such as wind, water level drop along the straits, 
water balance, topographic influences and the depth and width 
of the straits also must be considered. The scale of the 
motion and complexity of these factors make it difficult to 
determine the current system mathematically. 

Another main characteristic of these straits is the two- 


layer stratification of water which is associated with the 
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two-layer current system. This stratification is also a 
consequence of the water mass differences at the two ends 
of the straits. 

Detailed observations of temperature and salinity dis-— 
tributions indicate certain features of the diffusion of 
the Mediterranean water into the Black Sea. This water 
flows North through the Strait of Istanbul from the Sea of 
Marmara. It is shown by Bogdanova (1961, 1965) that this 
underflow of Mediterranean water into the Black Sea is 
present throughout the year and descends to significant 
depths in the Black Sea. 

This water exchange through the straits system fluctu- 
ates according to seasons, years and periods. This is very 
important for the hydrologic regimes of the Sea of Marmara 
and the Black Sea. Water level differences, given with the 
following expression, are mainly due to excess precipita- 


tion over evaporation and river runoff in the Black Sea. 


ne Zn — 2y (1.1) 


sea level differential between 
Black Sea and Mediterranean 


Ah 


Zp : water level of Black Sea 


Z 


M water level of Mediterranean 


Sea level differential (Ah) is always positive and fluctu- 


ates with time. It can be expressed in two components 


Ah = Ah + (Ah)' (le2)) 


a7 





where 


Ah : average sea level differential 
(Ah) ': fluctuating part of sea level 
differential [(Ah)' <« £(t)] 
According to Bogdanova (1965) Ah has a value 42 cm 
based on the long term observations. The fluctuating part 
of the sea level differential and total level differential 


are given in Table I. 


TABLE I. 


Monthly average fluctuation sea level differential 
and sea level differential. 


Month iit il vey ie VET vit 6X. UX.LUXE. UME 
(Ah) ' (cm) 1 0 3 (ape lS 6 0 “5-7 -7 —4 


Ah (cm) 43 42 45 49 53 57 48 42 37 35 soe 
(tBased on the data given by Bogdanova, 1965) 


The maximum change is 


2 


2a. ae 
xe (Ah) 5 Saye (Ah [ews = 9 cm/month, 


x 


observed during approximately June. This falls approximately 
in a period when continental runoff begins to decrease. The 
Fluctuating sea level differential increases from March to 
July. As a consequence during this period the upper current 


becomes stronger and the lower current becomes weaker. From 


oho 





August to November sea level fluctuation decreases and be— 
Bins to increase again from August to February. The upper 
current of the Strait of Istanbul becomes weaker and the 
lower current becomes stronger. 

According to Bogdanova's (1965) observations the maxi~- 
mum Ah value is observed during June and has a value 57 
cm; the minimum Ah value 35 cm is observed during October 
and November. 

The average sea level differential between the Black 
Sea and the Sea of Marmara is about 35 cm which is six 
times Moller's (1928) estimate (Gunnerson and Ozturgut, 
1974). 

There is a characteristic stratification with a sharp 
density transition layer over the whole area. The sharpness 
1S more pronounced at the North at the Strait of Istanbul. 

A cross section of sigma—t values for the Strait of Istanbul 
1s shown in Figure 3. Fresh warm water of the Black Sea 
lies in a surface layer 15-20 m thick at Istanbul. This 
water which is mainly river runoff and excess precipitation 
overrides the Mediterranean water which is cold and dense. 
The wedge form of the upper water shows clearly in this 
strait, but this is not true for the Strait of Canakkale. 

On the other hand, the wedge form of the lower water iS com 
paratively stronger and more pronounced in the Strait of 
Canakkale due to the downward sloping of the sea bed toward 
the North in the Strait of Istanbul. Average sigma~t values 


for surface and bottom water are 13.5 and 27.5. (Defant, 1961) 
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According to Gunnerson and Ozturgut (1974), the average 
salinity of Black Sea waters at a station 9 km from the 
Northern entrance of the Strait of Istanbul was approximate— 
ie 17.5 ° /00 for an observation period July 1966 to December 
1967. Lower values were observed (16-17 ° 7/00) POM 7 
to August 23, 1967. This lower surface salinity reflects 
maximum discharge of rivers into the Black Sea. Higher 
average surface salinities (17.5 — 19.0 ° 7/00) were observed 
for the rest of the year and extreme values (as high as 
Z5 °7/00) during winter months. These were due to the favor— 
able conditions for Mediterranean water to flow northward. | 
These conditions include small differences between the levels 
of the Black Sea and the Mediterranean and southerly or 
southwesterly winds over the Sea of Marmara and Strait of 
Istanbul. At oie midpoint the average salinity in the lower 
layer iS approximately 38.5 ° 7/00 and in the surface layer 
7.5 °7/oo. A salinity cross section is given in Figure 4. 

Warm and fresh Black Sea water can be identified easily 
from the temperature cross—section which is given in Figure 
5. Below this warm fresh water of the Black Sea cold inter= 
mediate water of Sea of Marmara can be observed. Mediter=— 
ranean water occupies the deep bottom layer of the Strait of 
Istanbul. 

The dynamics of the processes occurring in these two 
Straits is determined by the wind stress over the straits 
and the drop of the sea levels and densities at the end of 


the straits. The main driving forces are gradient of water 
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level, thermohaline forcing, wind force exerted on the sur- 
face and coriolis force. 

According to observations the boundary between the water 
masses generally does not coincide with the boundary between 
currents. The boundary between water masses has a greater 
Slope than the boundary between currents (Defant, 1961). 
This is probably due to diffusion processes occurring at the 
Current boundary and to the nature of the current (Figure 2). 

The average Black Sea inflow and outflow through the 
Straits of Istanbul according to the analysis of Moller (1928) 
are 6,100 and 12,600 m?/sec, respectively (Gunnerson and 
Ozturgut, 1974). Velocity is greatest at the sea surface 
and decreases rapidly with depth and laterally; it increases 
from North to South. Under normal conditions it is 40-50 
cm/sec at the entrance of the strait and 150 cm/sec at the 
South end. The lower current.is strongest in the central 
part of the lower water which is at about 16 m from the 
bottom in the Strait of Istanbul and 45 m in the Strait of 
Canakkale. These velocities are 100-150 cm/sec and 25 to 
10 cm/sec respectively (Defant, 1961). 

In general in the Sea of Marmara there is Black Sea 
water at the surface and Mediterranean water at the bottom. 
This is the normal state. Except for times when very strong 
southerly winds bring bottom water to the surface and cause 
Wind mixing these two main water types are both present with 
a very thin intermediate water layer. The Black Sea water 


extends to a depth of 15-20 m in the vicinity of the Strait 
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of Istanbul and becomes shallower farther south. There is 
intermediate water between these two water types which ex— 
tends at most to about 50-60 m. Below this the basin is 
filled with Mediterranean water as shown in Figures 6,7 and 
8. This water has a temperature 14.9°C, salinity 38.55 ° 7/00 
and Sigma-t 28.75, approximately in June. 

Surface water has a large annual variability. On the 
Other hand bottom water does not. Due to the described 
Oceanographic conditions bottom water is almost isolated from 
surface processes and has a large time-scale of variability 
compared to surface water. Average inflow and outflow rates 
at the straits give approximately 30 years for "turnover" 
time or "flushing" time for bottom water and 200 days for 
surface water. 

Based on the same rates a typical value for the vertical 
mean current |u| is 107 cm sec +. A characteristic magni— 
tude for horizontal velocity |u|] using the thermal wind 
relation is 4 cm sec +. Due to the large isolated reservoir 
of Mediterranean water of almost uniform characteristics in 
the lower part of the basin, the thermohaline circulation 


1s confined to the upper 50-60 m. 
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III. SELECTION OF THE NUMERICAL MODEL 


Regional circulation models may be separated into two 
categories. One type of regional circulation model is a 
diagnostic density model. This approach is based on an ob- 
served density field and wind stress distribution. These 
quantities, however, are not observed continuously and in 
an evenly distributed way. But based on average values of 
these quantities a velocity field is calculated diagnosti- 
cally. A diagnostic model, which seems at first easier, 
has drawbacks mainly of two types. First, it requires tak—- 
ing into account actual configuration of the basin geometry 
and, second, it requires sufficient representative observa~ 
tions. These two factors make this approach complicated 
and expensive. 

A second type of approach to study regional circulation 
models is to use a predictive model with an idealized topo- 
graphy and basin area. But this becomes more complex by the 
inclusion of more dynamical processes. 

At first glance it seems that diagnostic models are a 
better way of simulation of dynamical processes. But these 
observations on which such models depend are generally 
gathered by land based stations (e.g., wind data) and extra- 
polated to the area of concern. Furthermore, using average 


values brings up the question of how representative the 


22 





Besult is. Further, it 1s difficult to get a SeomtSPE SINE 
view of the dynamical process and long-term quantative 
analyses. 

On the other hand, although prognostic models are complex 
due to the dynamical processes involved, they are less expen=— 
Sive and can give a continuous view of yar epee which is 
generally based on an initial density field. At least as a 
first approximation it seems reasonable and economical to 
Study a regional circulation simulation. Present data are 
still used to estimate the initial density field and wind 
stress distribution in this kind of model. 

Due to the uneven distribution of density and wind stress 
observations it was mandatory to use the prognostic approach 
at the beginning. The procedure adopted is to use present 
data as efficiently as possible to estimate a realistic ini- 
tial density field and a constant wind stress distribution 


throughout the integration period. 
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IV. DESCRIPTION OF THE DYNAMIC MODEL 


A. BASIC EQUATIONS OF THE MODEL 

Basic assumptions which lead to important simplifications 
are: (1) the hydrostatic assumption, (2) the Boussinesg assump— 
tion with respect to density variations (density variation is 
neglected except where it appears as a coefficient of the 
gravitational constant), and (3) an implicit treatment of the 
mixing process by eddy diffusion. In addition, the non-linear 
terms are neglected in the momentum equation. 

With these, the sea is assumed incompressible, and density 
is replaced by a constant value everywhere except in the hydro- 
static equation. 


The governing equations based on these assumptions are: 


2 

eu OP 2 gu 
aE = a ee) oe enh LV (Get) 

Oo OZ 
wi Bey sayy 2Y~ tu (4.2) 
ot aoe M V 92° 
oP 
~ | (p Pg (4.3) 
gu OV OW _ 
ox t ty * 32 = 0 ee 
Go ig BE. aX) (4.5) 
ge 0H «P ap en ° 
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where the symbols have the following meanings: 


~ 


SS 


a 
= 


time 
horizontal coordinate meaSures positive eastward 
horizontal coordinate measures positive northward 


vertical coordinate measures positive upward with 
an origin at mean water level 


eastward velocity component 

northward velocity component 

vertical velocity component 

density 

reference constant density 

departure, from reference pressure which is 
pressure at reference density (= - 0,92) 
acceleration of gravity 

coriolis parameter (= 2 sino) 

angular speed of earth rotation 
latitude | 


coefficient of horizontal diffusion for momentum 
(constant) 


coefficient of vertical diffusion for momentum 
(constant) 


coefficient of horizontal diffusion for density 
(constant) 


coefficient of vertical diffusion for density 
(constant) 


0 () 


Ox 


0 () 3 () 


+v =~’ + w 


material derivative [= oe u ay ae 


introduces the effect of convective process when 
the stratification is unstable 


Si 





The sea is driven by wind and thermohaline forcing, these 
effects being introduced with surface boundary conditions in 
the momentum and density equations, respectively. The bound— 


ary conditions specified at the sea surface are 


iO 
K 35 0 
ou _ 
KY Oz 0 
y Zo) (4.6) 
OV oS 
AY Oz (OD 
Oo 
w= 0 





Heating or cooling at the sea surface is neglected over 
the period of integration. A linear north-south density gra- 
dient is specified in the uppermost level in the model and 
held constant during calculations. The x component of sur- 
face stress is taken to be zero and calculations are carried 
out with a uniform meridional wind stress. 

The boundary conditions specified at the flat horizontal 


bottom are 


JO _ 
K 35 = 0 
Z= —H (as) 
Uu=Vv=w= 0 


33 





The boundary conditions at the side walls of the basin 


are 
00 _ = — — 
Au Dy u V 0 VY 0, Ly 
(4.8) 
Orie eee = 
AL jy 7 U V 0 x 0, Lx 
where Ly : width of the basin ( 50 km) 
Lx : length of the basin(200 km) 


Boundary conditions at the straits are based on specified 
density and velocity fields which were taken from available 


data and annual average inflow and outflow rates respectively. 


Pee DESCRIPTION OF THE SOLUTION TECHNIQUES 

Since the model sea has a horizontal flat bottom and 
"rigid lid" approximation, vertical velocity vanishes at the 
surface and bottom. Following the method used by Bryan and 
Cox (1967~1968a~b), Haney (1974) and others, the ocean 
Surface was considered a balanced surface at which w = 0. 
In this approach the height of the sea surface is not obtained 
from the continuity equation; therefore, external gravity 
waves are eliminated. This also removes the vertical mean 
divergence from all scales of motion. Elimination of external 
gravity waves allows the use of a longer time step. In the 
case of long-term integration this is an important gain com 


pared to the loss of the prediction for height of the sea 
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surface. Filtering external gravity waves has a negligible 
effect on the density-driven part of the current. 

Due to the conditions on vertical velocity, integration 
of the continuity equation from bottom to surface requires 


the following condition: 


su. (OV 
poy 2 (4..9) 


where the meaning of the overbar is 


O 


Gi vertical mean [= = - () dz] 


—H 
and 


( )' indicates departure from vertical mean 
of any quantity. 
Also for some basic integral constraints this 1S a required 
relation to be consistent with boundary conditions. 
The "rigid lid" approximation, and (4.9), are accomplished 
by vertically integrating equations (4.1), (4.2) and subtract— 
ing the result from equations (4.1) and (4.2) respectively. 


Using present boundary conditions, the new sets of equations 





are 
ae 9p! Dee 97u! . 
oO dZ 
(4.10) 
2 t_7 
By 2 CP" ae DN a ee ' 
st poy MY Yt RY ToT 
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where 


u' eastward shear velocity component 
ee northward shear velocity component 
Pp! departure of P from its vertical average 
O 1 O O 
i= f (omegaan — ef ( f) (e-p,) gan) az (4.11) 
Z = Z 


The relation between shear current and total current for 


horizontal components is 


MeeVee) = ulkey,z,t) — u(x,y,t) 


(4.12) 
Vaan Zoe) = Vixsy,z2,C) — v(x,y,t) 
since 
lu} > > Jul and |v| >> |v| (45.13) 
equation (4.12) takes the form 
ah (Om ype) - GS Vig 2 ge) 
(4.14) 


iz 


I Vg ZC ) Maly 2, G) 


Thus total u, v component of horizontal velocity can be 
represented by vertical shear currents u', v'. In the 


following part primes of u' and v' are dropped. 
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The equations are non—dimensionalized by making the fol- 


lowing substitutions for new non-dimensional variables 


where 


(x,y) 
—> (Z) 


(u,v) 


mae 0) 


(P,P')= 


(t) 


(Tt) 


(o) 


(x,y) /L 
(z)/h 


L 
{PB , PT) (4.15) 
1 


scale velocity associated with density driven 


current 


scale depth, taken equal to a characteristic 
thermocline depth 


length scale 


wind stress [= MAX(t* (x,y), t” (x,y) ] 


Sigma~t [= (p~-1l) x 10°] 


surface Sigma~t at northern part of the basin 


surface sigma-t at southern part of the basin 


oF 





A scale velocity associated with the density driven cur- 
rent is defined by using geostrophic and hydrostatic scaling 
which were introduced by Bryan and Cox (1968a). When the geo— 
Sstrophic relation is differentiated with respect to z and the 


hydrostatic equation is used, the result is 


99 . ¢ 9U (4.16) 


Ne: 
a 


A scale velocity connected with the density distribution 


may be defined by 





= Gg Ap 
a > EL 2 (4.17) 
O 
where 
Ap north—south surface density differences 
(= Py 7 Pg) 


VV, °F (4.18) 


where 

wind stress 

L relative characteristic length scale 
for wind stress [= L/R] and Lisa 


length scale for wind stress and 
R radius of the earth 
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Substitution of these non-dimensional variables into 


(4.10), 


(4.11), (4.4) and (4.5) with some rearrangement, 


gives these equations in terms of non-dimensional variables 


of the form 


Pp! 





E E 2 
1 oP! H {2 V u 1 
Meee 2 vaty ri ty 
RO doy RO RO wy: RO LH, S 
1 
Oo et Oo O 
f (o-0 .) dn = a | ¢ i (o-o .) dn) dz (4.21) 
Z —H Z 
Oo 
f Ve V dn (4.22) 
Z 
E Z 
00 00 J0 i Zz K Vado 
Ua tV— twee e—Vor+r(ss) —_— 
ox oy OZ Pe K. Ro y2° 
6 (a) (4.23) 


where 


RO 


Rossby number which shows the relative importance of 
local time change term in the equation of motion with 
respect to that of the coriolis term Vy 
[= =] 
fL 


horizontal Ekman number which shows relative impor- 
tance of lateral diffusion and coriolis terms 
A 
_ M 
[= =] 
fL 


vertical Ekman number, gives the ratio of vertical 
diffusion of momentum to the coriolis term 


Ey 


{= — 


£h- 
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Pe Peclet number, shows importance of advection compared 
to diffusion of density 
Vit 
= a 


An 


V relative velocity, ratio of barotropic velocity to 


~ paroclinic velocity 

V 

D2 

[= 7! 

1 

HL. relative depth, ratio of basin depth to thermocline 
depth 

_ 4H 

[= Fl 


The model is governed by these six parameters assuming a 
value for L. is defined. All the parameters of the resultant 
run and constants of the model are given in table 2. 

Equations (4.19) — (4.23) are solved for variables u, v, 
P', wand o. The variables u, v, and o are predicted from 
(4.19), (4.20) and (4.23). P‘' and w are calculated diagnosti- 


cally from equations (4.21) and (4.22), respectively. 
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roo ae WD fe et 
Ky O 


A 


< 


< 
= 


Pope 


TABLE 2. 


3 
2 
2 


4.255 x 10. 
0.532 x 10. 
0.851 x 10. 
Or 0 

50.0 

0.1 

10 x i= cm 
Salk 107 cm 
8 minute 
iow (Gay 
40° 

27 day 
6.37 x 10° cm 
1000 cm sec - 
1.012 gm cm 
1 x 10° cm 

10 


1 


1.6 ene en 


B22 om* sec 


4 cm sec t 
a xX 10! cm 
8 2 —l 


0.5 x 10° cm™ sec 


0.4 x 10/ om? sec 


PARAMETERS 


AND CONSTANTS 


Rossby number 
Horizontal Ekman number 
Vertical Ekman number 
Péclet number ze 
Relative depth 

Relative velocity 

Zonal grid spacing 
Meridional grid spacing 
Time step 

Time scale 

Reference latitude 
Rotation rate of the earth 
Radius of the earth 
Gravity acceleration 
Reference density 

Depth of the sea 


Curl factor or relative length 
scale for wind stress 


Eddy diffusivity 
Eddy viscosity 


Characteristic baroclinic 
velocity 


Characteristic length scale 
Lateral eddy viscosity 


Lateral eddy diffusion coefficient 
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V. DESCRIPTION OF THE NUMERICAL MODEL 


A. SPACE AND TIME DIFFERENCING TECHNIQUES 

The space and time differencing schemes are similar to 
the schemes used by Haney (1974) for a general circulation 
model. 

In the integration of the primitive equations a staggered 
grid and centered space differencing scheme are utilized. 
The main concern in choosing the staggered grid was avoidance 
of a computational mode which appears when a centered space 
differencing scheme is used with an unstaggered grid. This 
computational mode in space is not present for a staggered 
grid. The staggered grid also has advantage in saving comput— 
ing time due to the presence of variables at alternate grid 
points. 

A leapfrog scheme is used for time integration. The 
starting scheme is a combination of forward and leapfrog 


scheme as shown below 


Leapfrog %At 


a 


n=*5 


n=0 n=l 


Forward At 


(5 time step forward) + (5 time step leapfrog) > Time step l. 
where 


n represents a time step 
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In most circulation models Matsuno scheme (Matsuno, 1966) 
is used as a starting scheme and to remove the solution 
separation which is present when a leapfrog scheme utilized 
exclusively for time integration. Since time integration in 
this model is not long enough to cause solution separation 
a cheme such as Matsuno is not used here. 

Each term in the momentum and Sigma-t equations is evalua- 
ted at a different time step to ensure linear computational 
Stability. 

Neglecting coefficients, the time differencing of the 


equation of motion corresponding to (4.20) is given by 


vett) 2 yl) 4 antipr(™) + pr lB); 
Sy eeriGins ez ut) | (5.1) 
where 
Ea pressure term 
ee friction term, which is composed of lateral 


diffusion of momentum and vertical eddy stress 
A trapezoidal implicit scheme is utilized for the coriolis 
term. Then, to satisfy consistency and linear computational 


Stability, coefficients Cl and C2 must satisfy the relation. 


CieteCZ =—1 


1 


The values Cl = 0.55 and C2 = 0.45 are commonly accepted 


for these constants which produce very slight damping of 
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inertial oscillations. Equation (5.1) and an analogous equa— 


(n+1) 


meon for u are solved simultaneously for the two unknowns 


ae), ae The time differencing form of the sigma—t 


equation has the form 


g(mtl) _ ,(n-l) (n—1), 


& 2at fapv '™) + ap Pl) 4 yp 


nt1, fee 


+ § (0 
Cc 
where 
ADV advection term 
HD horizontal diffusion term 


VD vertical diffusion term 


The leapfrog scheme is conditionally stable and the time 


step must satisfy the.relation (for two dimensional wave 


propagation) 
cee < (5.4) 
V2 
where 
C maximum phase speed of the waves present 


in the system. 


External gravity waves are removed and inertial oscilla- 
tions are rendered neutral by (5.1) and (5.2). Therefore 
only internal gravity waves are present to determine C. The 


phase speed of the internal gravity waves is 


Cc = vg"H (55:5) 
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Ag 


where 


qu is modified acceleration of gravity and 
given by the relation 
A 
g' = 4 (5.6) 
oO 
where 
g acceleration of gravity 
Ao a typical value for the vertical density 
differences 
O reference density 


O 


Assuming g = 1000 cm sec “, Apo = 10 x Og gm em - and 
a = eeol2 gm cm > gives a value for g' = 9.88 cm sec *. 
Integrations are carried out with an 8—minute time step for 


a grid spacing of 5 km, without any stability problem. 


B. THE FINITE DIFFERENCE EQUATIONS 

The numerical method is set down in terms of the non- 
dimensional variables. The finite difference scheme is based 
on a three—dimensional array of points with indices i,j,k. 
Time step is indicated by (n) as a superscript. In the nota— 
tion the letters i,j,k are always integers. 

Horizontal spacing is such that Ax = 10 km and Ay = 5 km, 
and uniform in the x and y directions. The horizontal grid 
pattern is shown in figure 9. Where horizontal velocity com 
ponents (u,v) are defined at circled points, vertical veloc 


ity (w) and sigma—t (0) are defined at dot points. 
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J:1 
=1 





eo u,v Point 
e T,w 
Point 


Figure 9. Placement of variables on 
horizontal grid plane. The open cir- 
cles denote the definition points of 
the horizontal velocity components 
(u,v) and dots denote definition points 
of (o,w). The distance between adja~ 
cent grid point is Ax = 10 km and 

Ay = 5 km in x and y direction, res- 
pectively. IM and JM have values 21 
and ll respectively. 
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For an arbitrary interior point, where (u,v) are stored 


longitude (x) and latitude (y) are defined by 


x +5= iGo) J J = Mp2. es op M1 
(5.7) 
y, +52 H+ Gv oe epee 


The vertical pattern of variables is shown in figure 10. 
The vertical index indicated with k and the vertical velocity 
w are given along the surface and bottom and along the layer 
boundaries. The variables (u,v,o) are located at the kth 


level. Depth of each level is defined as 


y, =_ — aan 
1 2 
(Src) 
m=k 
Zh = Za) = : (AZ ro AZ -1)/2 Ke =" 22 es wl 
m= 


Layer thicknesses for each level are given in table 3. 


Table 3. LAYER THICKNESSES 


= 0.25 0.25 0.3125 0.6875 1.0 1.5 6.0 15.0 25.0 <& 


To resolve the surface layer more accurately, four levels 
are located within 20 meters. Close to the bottom a coarser 
grid is used. Resolution in the surface layer is greater than 


in the bottom layer. 
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1 U,V, ¢ 9.5 
ay Eg 

2 Bava 7.5 
eee ee eee ee 

3 U.V.¢ 1 2.5 


|# 








4 UV eee lO 
oe 
5 U.vie 40.0 
er eae 
Cee ee eNO LO 
bly w 
K- 7 U,v,T 10 0.0 
7 eee WO 
ee Re C=*=BBFD':1 
ee eon ee ee 
KM=9 Uw 700.0 
9%, . oe sd ee OOOO 


Figure 10. Vertical structure of the 
model (u,v,o) are defined at integer K 
values and w is defined along the layer 
boundaries. 


48 





The finite difference form of the hydrostatic equation 


has the form 


P... = (o-0_) AZ 


“Tapl Oi J, 3 4s 


+ (o-0 ,) AZ) 4 k=2,3,KM (569) 


ijk ~ Paij,k-1 ijk 
and the vertical average of the departure pressure- is 
calculated with the relation 


k=KM 


Pi; = 2 P isk AZ (5.10) 


The departure P' is calculated from 


Pe, = P (52.12) 


i5 ijk ~ Pa3 
fieen 1S the finite difference analog of (4.23). 
The finite difference approximation of the continuity 


equation, in which the vertical velocity between the levels 


is calculated, has the form 


Ss eik+h = ee + [U, + 44k AZ, (5...82)) 
where 
— os 
Smee ee Ox jee * Ux jk ik 
(5:43) 
+vs. ’ 
a ijk = Bor it y ee 
and 
u cee = + (y — u ) 
x 1j+%+k Ax ~~ 1i+%5 1—4s" j+3sk 
(5.14) 
V = iy mn) ) 
y it+sjk Ay * j+%s jms it%sk 
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which are analogs of ($2) and (su) respectively. Gradients 
of horizontal velocity components uy and vy are defined at half- 
integer grid points west and south of the u and v storage 
points shown in figure lla. 
Integration from the surface downward, uSing surface 


boundary condition w=0, is done by 


k 


Wi skeh = » (U,, as Co AZ (>.>) 


The formula for the finite difference approximation of 


the equation (4.19) by which u is predicted is 


n+1 n—1 Pp! — Pp! 


(Se = — 1 ,;_iti jtl i jt 
Ze 1+%j+%k 2RO Ax 
t oe 4 _ 
» Ath GT Phin y SH (Px ind jats 7 Oe bins 
Ax k RO Ay 
» dy ds j4a 7 Sy iesjynea , Bv (Be koe 7 Ye kt nol 
Ay k Ro AZ, it% j+% 
iL n+1 n—1 
a. oe i 
a +14 +3, 
ee he hy Dean (5.17) 
u —u 
k=] k 
and eee ey 
Maks AZ iti +s (5.18) 


Equations (5.16) and (5.17) represent (5) and (se) and 
are located half grid distance south and above the velocity 


(u) storage points, respectively, as shown in figure 11 a-b. 


5:0 





2° S3uTod obezojs oO oy} saAaoge 
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aaoqe sjutod ptib JTey ye oze “A pue “n *o pue A ‘nn JO sqUaTperb 


BU SUTJAP YOTYM SETqeTAeA FO SUOTREOOT TeROTAIADA °q 


*peuTsep ere soTqetTAeA ssoyy 
JO sqZueTpeIb squtod x 3W °O FO suOoTzeDOT BUSSeAdeA sjOp ‘“SaeTOATSO 
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The formula for the finite difference approximation of 


the equation (4.20) in which v is predicted has the form 





— 0 paar t 
aa eee itl je Party 
2At ith jtsk 2RO Ay 
0 — 
Mm ij+1 Piijn A By vx itl jt — Vx ij+s 
Ay k RO | Ax 
P Vy itkjtl — ‘y i+%j) (n-1) 
Ay k 

E,, ie ait Viet 

RO AZ, i+45j) +4 

ve I, 

pr sr WY _ +t n+1 n—l 
* Ro H * it RO [Clu + C2u" Vyasa (5.19) 

where 
Ve 
i+1ls its 
Vx i+l Ax ) 4 sak (5.20) 
V ae 
a k 
5 aes AG ) ees (5.21) 


Equations (5.20) and (5.21) are analogues of (Ss) and 
(=) and are stored south and above the velocity storage 


points, as shown in figure lla-b. 


SZ 


The finite difference approximation of the sigma-t equa- 


tion has the following form, 


ert a og i, 2 Ox ite Ox ej 4 
2 Nt ijk 13k Pe MG 


Oy itkjth Cy itkj—s,n-1 
+ Hate ne 


Ay 
i Kid1+, (2 kt ~ C2 k~g) n—1 
K, Ro “v AZ ") 44355 +35 
n+l 
+ 6 (0 isk (S.202) 
where 
= ——((u Pave = i= (6 5.4.) 
fir) K 4Ax it%j+% 1+%j— ij i+lj‘k 
e (Uj 454% o Uj_35-4) (F545 ree 
+ ar + Vv )-(o Ce) 
4Ay itsjt+ Deon) aes I5—1 Lak 
Beare Viegas) FG 5a 7 945" x! 
eee Cw. )e(o 1 ) 
2AZ, 1jk—5 aK ijk-1 
~ (Wi 543s) ag F Fi5Ker)! Pay 


5 





and 


pee eae al Oa’ Gk 
Oo = — (o — Oo ) 
y ijt+%sk Ay j+1 Lk P5244) 
OF iik+s ~ = CP me ee 4 
J 2 k+45 Zz J 


As shown in figure (ll a-—b) Os and o and o 
2 


y j7% Zk—%5 
are defined at half grid distances south, east and above the 
Sigma—t storage points. 

. Before proceeding to a new time step convective adjustment 
is applied. It is difficult to state the convective adjust— 


ment mechanism. With this mechanism the density structure is 


forced to remain stable. It may be expressed by 


é = (-) (Go25) 


In case of unstable stratification, vertical mixing be- 
comes effectively infinite. At each time step this infinite 
mixing is included in numerical solution by testing Sigma-t 
meerrle for unstable lapse rate (Op > O%)- 


tion exists new values of sigma-t for those two layers are 


If this condi-— 
Set equal to the vertical average sigma-t of the two layers 


instantaneously. This process is repeated until complete sta- 


bility is reached for the entire layer. 
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C. SPECIFICATION OF BOUNDATION CONDITIONS 

As shown in Figure 3-1 06,w are variables which are de- 
fined at lateral boundaries; also horizontal gradients of 
“{u,v) are defined at these lateral boundaries. The boundary 
conditions are satisfied by the following relations which 
are given for eastern and western boundaries. Similar con— 
ditions exist for north and south boundaries with the excep- 
tion of open boundaries. 


Normal components are set equal to zero by 


ee _ 2U)- 
X in] Ax’ i1=4 

(5.26) 
. a 2u) —_— 


Se 


In the momentum equation zero—slip condition is applied 
to the valocity tangent to the boundary. On the other hand, 
when computing the advection term in the density equation 
free slip condition exists for the velocity tangent to the 


boundary. This implies the following conditions 


1=1*s 
(5.27) 
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Definition of vertical gradients of (u,v) allows the 
writing of the vertical stress term in flux form easily. 
Boundary conditions at the surface and bottom are defined 


in the following manner: 
The 
uu. k=5 = 0 (5220) 


Zonal wind stress is assumed zero and only a constant 
meridional wind stress is defined. 


V_L 


eee 
Ve k = 5 EL ive (5.29) 


On the other hand it is assumed that horizontal velocities 


vanish at the bottom. The following conditions force veloc— 


ities (u,v) to reach zero at the bottom 


- = 2u 
Z k=KM+4;  2Z!k=KM 


(5.30) 
o ee 
Z k=KM+3; AZ! k=KM 


The boundary conditions of zero mass fluxes across the 
Side walls and bottom are satisfied by imposing the following 


conditions: 


(ou) gaa, = ~(0U) 23721 5x 
(5.31) 


(OW) sere ~ 80) sory! ay 
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Similar conditions exist at north and south boundaries. 


In the diffusion term zero flux and insulation are 
satisfied by the following conditions: 
0 ———  O 
*i =k *i=3/2 
(5.32) 


Ox = “Oy 
1=IM+3; i=IM—s 


where 05. is defined at half grid distances west of the 
sigma—t storage points, as shown in figure lla. Similar 
conditions exist at north and south for meridional gradients 
of sigma-t, with the exception of open boundaries where 
Sigma—t is defined outside the domain. 


At the sea surface downward flux is absent. 


eo ke = = 0 ash 


The sigma-t equation is solved for only the upper six 
levels. It is assumed that sigma—t for the rest of the depth 
range has a value Our which is held constant during integra— 
tion. This enters the system in the calculation of the 
vertical gradient of sigma-t one half grid distance below 


the sixth level, 


Oo = 


(5.34) 
ke =6% 


and O. has a value 28.5 in sigma—t units. 
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VI. RESULTS 


Computations are carried out for different values of the 
parameters with an integration period of 24 hours. Parameters 
of these runs are given in Table 2. Results of the model are 
given for the fields of horizontal velocity components (u,v), 
vertical velocity (w) and density (sigma=—t) in dimensional 
form. 

The horizontal velocity field at 2 = —2.5 m is given in 
Figure 12. The horizontal stress exerted on the sea surface 
by southerly wind causes the surface water to move ina 
generally North-east direction. This is a consequence of 
movement of surface water as an Ekman layer, drifting to the 
right of the wind stress in the Northern Hemisphere. Since 
the curl of the wind is zero divergence of the Ekman drift is 
also zero. But horizontal velocity is strongly convergent and 
divergent near the boundaries. Northeastward drift takes 
place over the entire basin at this depth. The strongest flow 
has a magnitude 20 cm/sec in the central part of the basin and 
diminishes toward the boundaries. Near open boundaries 
velocities are smaller compared to the velocities at neighbor- 
ing grid points. The southward flow at these points repre- 
sents the exchange with the Black Sea water at the north and 
Mediterranean water at the south. The magnitudes of these 
currents are 5 cm/sec and 2 cm/sec at grid points just near 
the entrances of the Straits of Canakkale and the Strait of 


Istanbul. The applied meridional wind stress has a value of 
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0.3 dyn on” which 1s not strong enough to move surface water 
to the north as it does at neighboring grid points. 

In the second layer z = ~7.5 m the horizontal velocity 
changes direction to the right and becomes weaker. Horizontal 
velocities at that level are shown in figure ae: they have a 
magnitude 5 cm sect. Flow adjacent 'to the straits is still 
southerly eee a magnitude 5 cm sec + and 3 cm sec at 
the vicinities of the Straits of Istanbul and Canakkale 
respectively. 

The circulation pattern at lower levels differs from the 
upper level flows both in magnitude and direction. These 
differences can be observed in the circulation pattern at 
the z = —40.0 m level (Figure 14). The pattern is not irregu- 
lar compared to the upper layers. 

The northeastward surface drift in the first layers is 
accompanied by a westward and partly southwestward flow in 
the bottom layers. A northward flow adjacent to the Strait 
of Canakkale brings Mediterranean water into the basin. A 
Similar northward. flow of Mediterranean water does not appear 
at the northern strait due to the shallow sill depth at that 
point. 

Flow in the layers close to the bottom becomes weaker. 

The influence of the circulation at the straits on the general 
pattern is small. The circulation pattern at a depth 700.0 m 
is shown in figure 15. The maximum magnitude is 0.033 cm sec +. 

One advantage of the numerical modeling study is that 


it gives estimates of important oceanographic variables that 


60 





d °(ToejZouU S°L 
‘TaAeT puodes TOF SAOAVDSA AQZTOOTSA Te_QUOZTAOH “ET SANHTA , 








ee gee eh eet ee et! ty! Ve Pe . 
Sy AT RDO ZR RR tS AG 
Ny AZ SZ OE ep nh EEN 
Sg RTE Sk ee gp eh gg 
EB OR Ng Rg ed eg oe OE 
EE hg Ng Kg gp gh gd a 





ALP Ff AVN LN f-AAAAaAH SAKA 


61 


*(z930u 0°07) 
‘TOAST UAFTJ AOF SAOQZQODA AQTOOTSA TeRUOZTAOH “PT sanhbty 


ee te er Ye re a A AT 


T 


0 


(eee pent, Oo, Oe 0s GO Oe Oe OL, Bs p= f l 
a cence a =e Cy Cece ee Gees mE geet pe i ah 
ew oo eo Oe aaa aaaiaiadgoatdaoa ggg A 
op o-~ o—o-H2_A@a agai a giogiaiod og 9g g- g- 
gegen ee ee Ce Doe Can Dee Coe DD OL fg yg 
a or pe ll hoo 


wee eee Loe ooo eo be oe Oe g o 


We oe i of ald Lod 





62 





*(aZaz0W 0°O0L) 
‘TaASeT YRUTU AOJF SAIORZDSA ARTOOTSOA TeRUOZTAOH “*GT eanbty 









Ue a ee ee ae a ae Se J 
APP AAAAAAAA PTS YAP PPA 
Af PAAAAAAAAS AYaAASA AA 
AA PA AAA AAA AA AYA AAA AA 
po A tet 6 8 8 88 ee eg 5 et 4A a 
APA AAA AAA AAA AAA AA PA 
AA AAA AAA AAA AAA BAA APA 
AAP PAA AA AAA AAA BA A A A Pf 4 
eee) ae AA A A AA AP Ff 
MILLIS PLEAS Bis 1 lay 


63 





Cannot be obtained by measurements. One such variable is 
the vertical velocity. The vertical velocity at ene base 
of the first layer is given in Figure 16. The magnitude of 
the vertical velocity differs markedly between the boundaries 
and the interior. Rising motion (upwelling) is taking place 
in the south and southwest and sinking motion (downwelling) 
1s taking place in the north—northeast. Vertical velocity 
is mainly determined by the Ekman drift current intersecting 
the boundaries. Southerly wind stress produces upwelling 
and downwelling at the south and north, respectively. Up=- 
Welling extends to the north on the western side and down 
welling to the south on the eastern side. The maximum magni- 
tude of the vertical motion is 358 cm day -. Strong upwell= 
ing and downwelling take place adjacent to the straits. This 
is consistent with the nature of the horizontal current 
system in the vicinity of the straits. 

Vertical velocity at the bottom of the second layer, 
Z = -—7.5 m, is shown in Figure 17. The general pattern is 
Similar to the vertical velocity at the base of the first 
layer. Strong vertical velocity is present along the lateral 
walls to the southwest and northeast where it has a magnitude 
+ 300 cm day +. With depth the effect of wind on vertical 
velocity decreases due to the vertical stability. The verti- 
cal velocity pattern at the base of the eighth level is 
Shown in Figure 18. 

The sigma=t pattern at the first level, 2.5 m is shown 


in Figure 19. This pattern at z = —2.5 mis specified by 
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the surface boundary condition. At this level the surface 
sigma—-t is a function of latitude only and held constant 
during the eoretien period. 

Deviations from the surface pattern begin to appear just 
below the first level at 7.5m. The Sigma-t pattern for the 
second level is shown in Figure 20. A close relationship 
is present between this sigma—t pattern and the vertical 
velocity field shown in Figure 17. At the north and northeast 
Sinking motion brings in low density and water and at the 
south rising motion brings in high density water to this 
depth. This pattern iS more strongly tied to the surface 
boundary conditions on sigma-t than at other levels. There 
is a low density water pool at the north. As a consequence 
of upwelling a high dense water pool is present at the south 
extending in an east~west direction. Magnitudes of sigma~t 
are 16.5 and 18.0 at north and south respectively. 

The Sigma—t pattern at the fourth level z = —20 mis 
Shown in Figure 21. The eRe CeCe vertical motion appears 
as a low density water pool at the northeast. Below this 
level z = —40 (not shown), Sigma-t is uniform and has a 
value 28.5. This is consistent with the sigma-t pattern 
given in Figure 6. 

Further insight into these results may be gained by examin- 
ing the meridional and east-west cross sections of sigma—t 
and the (u,v) components of the horizontal velocity. There 
is also the possibility of comparing the predicted sigma—t 


pattern to the observations. 
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Meridional cross sections of sigma-t at different longi- 
tudes are shown in Figures 22-24. The general patterns are 
Similar. High density water rises to the surface and strati- 
fication becomes shallower further south. Below 40 m density 
is constant. This is consistent with the observed sigma-t 
field. . At the west and south very strong upwelling and weak 
downwelling at northare shown in Figure 22. In the central 
region there is upwelling in the south and downwelling in tke 
north as shown in Figure 23. Further east the effect of 
strong downwelling is shown by ene rate of tilting of the 
1sopycnals at the north compared to the south in Figure 24. 

Meridional cross sections of zonal velocity in the east, 
central and western parts of the basin are shown in Figures 
25-27. The zonal flow at the surface is eastward and changes 
direction with depth. Zonal flow is westward at greater 
depths. Strong flow is present at the surface and decreases 
rapidly with depth. Variation with longitude can be observed 
by comparing the figures. In the central part zonal velocity 
has a magnitude 16.0 cm/sec at the surface and decreases 
both laterally and vertically. Below 20 meters flow is west- 
ward and has a magnitude 0.6 cm sec *. 

East—west cross sections of Sigma-t patterns for three 
latitudes are shown in Figures 28-30. A regular stratifica-— 
tion is deformed by upwelling and downwelling of the eastern 
and western sides of the basin (Figure 29). Departures are 
present in the north and south, as shown in Figures 28 and 29 
respectively. Intrusion of less dense Black Sea water can 


be identified in Figure 30. 
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The vertical cross sections of the meridional velocity 
which correspond to the previous Sigma~-t Sections, are given 
in Figures 31-33. There is a northward flow at the surface 
in each section. The cross section near the southern bound— 
ary shows a southerly flow below the surface layer. Flow 
changes direction at about 20 meters. In the vicintiy of 
the Strait of Istanbul the effect of the S6pen boundary is 
indicated by a weak northerly flow. There is a southerly 
flow below that and at depth about 30 m another reversal 
takes place. Weak northward flow is associated with the 
southerly wind. The southward and northward flows at depth 


are consistent with observations. 
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VII. CONCLUSIONS 


In order to investigate the dynamics of the Sea of Mar- 
mara a.circulation model is formulated based on the hydro- 
static and BousSsinesq approximations. Since very little is 
known about the turbulence characteristics of the ocean or 
adjacent seas, viscosity and conductivity are replaced by 

new terms representing the contributions of the smaller scale 
motions to the exchange of momentum and density. Horizontal 
eddy transfer of momentum is accomplished using a constant 
eddy viscosity and diffusivity coefficient. For the exchange 
of momentum and density in the vertical a constant verecenl 
diffusion coefficient also is used. The vertical mixing of 
density is enhanced by the use of an instantaneous adjustment 
mechanism to eliminate unstable lapse rates. A lateral mix- 
ing coefficient for density is less than that used for 
momentum. These coefficients should be eae een ie small 

ee as not to obscure lateral transfer of the quantities, and 
a non zero momentum coefficient is needed to satisfy the 

zero Slip boundary condition. In this study constant values 


of the coefficients are 


~or xs OL em*sec (density) 


*a 


and 


An 5 x jen sccm (momentum) 
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The vertical diffusion coefficient which is difficult to 
evaluate either by measurement or in principle, is assigned 
1.6 and 3.2 cm?sec > for density and momentum respectively. 
Other parameters and constants are given in Table II. 

Further study of thermocline penetration, water exchange, 
Pertscal mixing processes and current gradients might provide 
a better possibility for matching these coefficients. 

Although many important factors, such as the irregular 
bottom topography, bottom friction and non-linear terms in 
the equation of motion have been omitted, the results from 
the model are generally encouraging. More realistic results 
may be achieved by incorporating these effects systematically 
in further studies. 
> Integration in the model is carried out with a constant 
southerly wind stress which was chosen arbitrarily. Although 
the scale of the motion is small compared to a characteristic 
length scale for meteorological features, wind stress may be 
Significantly variable over the sea due to the geographic 
location of the sea. Also seasonal variations of the wind 
may have a conSiderable effect on the vertical mean part of 
the current. which is neglected during this study. 

The model does not separately predict temperature and 
salinity which are sometimes important quantities. Diagnostic 
calculation of the density with the aid of predicted tempera-— 
ture and salinity fields can be done with a very simple 
modification of the model. This would require a calculation 


of the heat and salt fluxes at the surface and meridional 
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salt and heat separately in the model. It would be more 
difficult, however, to define boundary conditions for these 
quantities at the open boundaries. It is also difficult to 
Simulate seasonal variations of these quantities at the 

Open boundaries. Once these are defined properly, based 

on accurate observations, there are no inherent difficulties 
in simulating the dynamical processes in the sea using the 
general principles involved in this model. Even though this 
model depends on the particular hypothesis used for hori- 
zontal and vertical eddy EeaneeOae of heat, salt and momentum, 


it is clear that valuable studies on small scale water bodies 


can be made. 
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APPENDIX A 


COMPUTER PROGRAM DESCRIPTION 


The computer program iS written in FORTRAN IV and was 
used with the IBM 360/67 computer system at the W. R. Church 
Computer Center, Naval Postgraduate School.’ 

The overall program is divided into two basic subpro= 
grams (1) the main program and associated subroutines (2) an 
access program which draws and writes the results of the 
first program. 

The main program consists of nine subroutines which 
calculate different terms of the equation of motion and the . 
density equation, pressure, vertical velocity and changes 
variables for next time step. 

FORTRAN IV symbols for the primary program and a brief 


description of the subroutines are given below: 


UMI Zonal velocity component at n=l time step 

U zonal velocity component at n time step 

UAL Zonal velocity component at n+l time step 

VM1 Meridional velocity component at n—l time step 
U Meridional velocity component at n time step 
UAL Meridional velocity component at n+l time step 
SGM1 Sigma-t at n—l time step 

SGMT Sigma~t at n time step 

SGMT Sigma—-t at n+l time step 
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ADV 
ADV 


ADV 


PBAR 


UXG 


UZG 


VXG 


VZG 


SGX 


SGY 


SGZ 


TX 


vy 


DZ 


DZ 


BB 


BO 


CD 


DX 


IY 


RO 


Total advection term in density equation 

Local rate of change of sigma-t 

Pressure 

Surface pressure 

Vertical average pressure 

Vertical velocity 

Gradient of zonal velocity in x direction 
Gradient of zonal velocity in y direction 
Gradient of zonal velocity in vertical 
Gradient of meridional velocity in x direction 
Gradient of meridional velocity in y direction 
ee aai cnt of meridional velocity in vertical 
Gradient of sigma-t in x direction 

Gradient of sigma-t in y direction 

Gradient of sigma-t in vertical 

Meridional wind stress | 

Zonal wind stress 

Layer depth 

Layer 


Depth dependent velocity at the Northern 
fictitious boundary 


Depth dependent velocity at the Southern 
fictitious boundary 

fictitious boundary 

Sligma—t defined outside the domain at North 
Sigma-t defined outside the domain at South 
Horizontal grid spacing in x direction 


Horizontal grid spacing in y direction 


Rossby number 
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EH 
EV 


PE 


HH 

DT 
SUBROUTINE 
SUBROUTINE 


SUBROUTINE 


SUBROUTINE 


SUBROUTINE 


SUBROUTINE 


SUBROUTINE 


SUBROUTINE 


Horizontal Ekman number 


Vertical Ekman number 


Péclet number 


Relative velocity 


Vertical eddy diffusivity 


Vertical eddy viscosity 


Curl factor 


Depth of the basin 


Time step 


PRES 
VERW 


ADVEC 


HDGRD 


Die © 


SIGEQ 


HVGRD 


HORV 


Calculates pressure 
Calculates vertical velocity 


Calculates advection term in the 
density equation 


Calculates gradients of the sigma~-t 


Calculates local time rate of the 
change Sigma~t by subtracting advec— 


tion term from the diffusion term 


Calculates new sigma~t and makes 
convective adjustment for unstable 
lapse rates. 


Calculates gradients of the hori-~ 
zontal velocities 


Calculates horizontal velocity 
components, u,v. 


A descriptive flow diagram of the program is shown in 


Figure 34. 
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Figure 34. Descriptive flow diagram of the program. 
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Figure 34, Continued 





35 





Yes 


! i 
BB, CA, DZ, TZ, TY, U, V, UXG, UYG, VXG, VYG, UZG, VZG 


{ 
t 
No | 


elelety ye UM1, VM1, UXG, 
Yes 


HV 
gv 
HORV 


D4T,DZ,ADV, UXG,VXG,UYG,VYG, UZG,VZG,TX,TY, U,V,UA1,VAL 


. 
t 


i: 
No |! 





HORV 






No 


D2T,...,UML,VM1,UA1,VA1 





Figure 34. Continued 


96 





= | 


DT, .ee,UM1,VM1,UA1,VAI1 


gv v 
HVGRD 
eoe0e ,UA1,VA1,USB,UYG. ee 
DZ, UXG, VYG,W 
No ‘ 


Yes 


UA1, U, VA1, V, SGP1, SGMT, SGP1, UMl1, U, VMl, V, SGMl, SGMT 






U, UAL, V, VA1, SGMT, SGP1, UMl, U, VMl, V, SGMl, SGMT 


Figure 34. Continued 


97 





YES 


Figure 34. Continued 


98 





uUija 
>nMQ 
wwe) 4 LL J 8 
UW of of 
AAA <_« 


ome 
ZN wd 
rng <I pong 
wv = f~LL 
aqa7a<at 
zOoa@a 
Oe Nike 
dome > 04 (/) 
<a 
Zz, ad ©) a 
> LD < 
=O 2 
=m x Co 
Om us 
O<te 
Ww = 
ULy o4/) 


IlerwO 
~-<rr 


<Irnuw 
<I 
MZAO 


LL. =e fa 
TO TzZ 


CL <T UL ms 

Owuoa 

On<dar< 

Zc tw 

Owiard 
x 


MNeMwi 


Tiki >5 
-Oawn 


WNOOO 


IelO eae 

D> enn 

@ ew 0 OO 

am eee mS & —_ 

DOW aeQi tt Oo 

eel | ld ~— (/) 

OOnW! eae <w 

eo NAIC —~ el 

NANO & O09 = LO) 3 

wee we /) CD = <~") 
a> eYD>> Cl) = 
re om © & & mO TL 


OND RONDO tem 
rte ord ee eo em ww OY 
NOOaAOOQO wll 
wate omtaeet ee je > 
OL em we © em & Med HM) <{LOmrAl 
LOON Nett —m OOS & 
QNASKNAIAIC’G FOO— @ eff) 
OL wee te ged wee ee ee ee we ON LO) LO fh @ 
AS XO DA ww TILILOO 
ZeZOOXKMOAWNXODESIL 
OD>2NM DN > ASO KO MAJIK 
= nee weet ene Sea? 
N22 ZZ2Z2Z222 22 Re 
ZOQQOOQVOOO0O0O0OONCST I 
WEL SSAA LT TTI STS 
2S LASALLE SASS Taare 
—ONdQOQVI0 0ON0QVWONI0OO 
OVO QOOQ VO OU UOLLLLLL EL 


MINOW 
MAIN 


X, 1F8.3) 
CALCULATION OF THE LAYER THICKNESSES 


READ TRE LAYER DEPTH 


WO 


p1=2148,2) 


THIS DO LOOP CALCULATES INTERMEDIATE DZ 
READ MERIDIONAL WIND STRESS TX 


C 
C 
C 
C 


—O 
r= om 


—) 
1 | 
=~ = 
oo @ 
SS Sh 
mf ome 
a wee? 


<>< 
hb 


a aes? 


AE ey 
uvu 
a= 

a & 


yu 


ag? Cag 


aA<q 
eat S) 
ew 





PARAMETERS OF THE MODEL 


C 


GJ 

iN 

tL | 

OvU 

a= pet 

CECT LLY end vd uy 

O eNO @ @ OHIO @ 
em «© COCO « eO 
V#8NOoo i =—oOown 
Ron nth nh 
MM oh Ld Sd KO 
—I<IWWwa>o> cool 


HORIZANTAL VELOCITY FIELDySTARTING AT REST 


C 


POINTS 


oo 
DWHNoOo 
one @ 6 
set OS 
ntanw i 
YD ot ore LL LL La 
we D>) 
HA © aS Za 
OOO 2 7a 
(NNIA @ @f--}-- 
ond 
QogOo~~""D 00 
OOaDd>O0OUN 


VELOCITIES AT THE FICTITIOUS GRID 


CT) No emt 


ooo 
NAN 


C 


100 


iy iy 
ee ee ff 
FN @¢@ ee eMieted @ 6 6 
| 1Onmgit 1ooo 
eo etka 
Ry LOT ay, <i, GOO, LR AE, LO TR, Ry 
KANO PMO a CUO PLUVO 
fog? Soge? Seg? aus? Sage? Meee? age? Sone? “qe? “age? ee 
QwMawaaaataada 
OoOMODDIGIIVIONO 


D SURFACE PRESSURE 





OO 


uN 
N 
=~ 
@©ee9e#e? @ 
Stet PPhy hr 
SNH HW WW 
S Ce 2 6 SS ee ee 
eee NI SPUN 0 
Hnhewwnweenen a 
DD DOLL LS 
ene eo oe a) 
ND [= 04 04 ng mt Od Od SS” 
QO wees we wee we wes wes Pd Od 
alee is =e Rie 
LARALAALS 
OOQVOOOOOU0 
QAQIMNNMNWNMNOO 


r=0 
Oo 
= 


-0O. 1*(J-1) 


LY LVL 
N NN 
om uN mo 


eesee@ee0e0e288eeeeft @ @ 
CO) PLL SO Pe SPUN OO R= Pe Poh 
| | | | | | | 
iil Tig, CP, PMR IT, IAB Ey GM) A py ay GR 
PION PUYN HHO SPU 
ea ee egy ee? MP aay Cael Sa? tag MP aay eg? 
WOOQQOUOWOOOWdId 
DOAMDMOADQIVYUO VU) 


ORE 


FIELD 


INITIAL VERTICAL VELOCITY 


=m 

Dal 
aeren e@ 
ste Hto 
wu 
~~ a= 
ww 
Ooo « 
LA LALA => 
NQor oa 
= 
WwoOoOdO~ 
QOO= 


© 
uty 


Oo 


= DT 


TIME STEP 


C 


le 92 E-4 


DT= 


jl i NG INTERVAL OF THE RESULTS 


WRITI 
K1=90 





(TdOS*TWOSS AQVSZ0%1L20)0Z9IS WWI 29 
O0c6 OL 09 
(TdOSSLWOSSAGVS 70%L70)039IS 11V9 19 
Oc6 OL O09 
(TdOSSTWOSSADVS 70% 10)039IS 119 
c9 OL OD (2°03°NjdI 
T9 OL 09 (T°BAPN)II 
SdLVeY ASdV1I AIDGVLISNN HOS 3211508d AVIILYSA LSNFAV 
QNV IW9S MAN S3LVINIIVI O3S9NIS ANILNOYINS 
(AQVSZ9S* ANS £X9S*70)033I0 T1V9 BT6 
4340339 Q31L1V INDTV)I 
SI LVHI Wu3L NOILISAQV SHL OL SQQV ANV WHsl 
NOTISNASIO S3ALVINIWI O4SSIGO SNILNOYITVS 
(AQVS IWOS fVOS aatMSAfn§ ees eas. 5S 
(AQVS IW9SfV9SAESMSASNSGD SOE 70) 93AGV 1199 
6S OL 09 (2°0A°N) AT LT6 
NOILVNOS ALISN3Q 3HL NI Weal NOILISAQV SALVINIIWI JSACV ANILNOYNS 
ee eee ee ee eee eS BS 
(Z9SSADS*X9OS § TWOS £V9'G9£09*0E*Z0)J0NNGH 11V9 
8S OL O09 (T°DAPN)AI TS 
SLNSIOVYS ALISN3SO WOIILYSA ONV IWINVZIYOH 
SSLVIANDITIVI GQYSQH ANILNOYINS 
NOTLVNOZS ALISN3Q YwOd NOTLNIOS 
N (8)3LI1YUM 99 
TS OL O09 
279704 OD Ee ecces Ket 
99 OL QD (eH°d3°N)ITI 
99 Gl O9 (THe O3SN)II 
XH’ T=N ¢S 00 
ALTIOIZA TVIILYSA ONV AYNSS3Yd BJHL AO 
NOILVINDIVI GNV LWOSSASN JO NOILIIGIYd 
O8T=X > 
GET=2y 


OVO 


102 


OO 





CULATES PRESSURE 
UR 


LCUL 
SURE 


In 


»ADV) 
AySGM1,0Z,ADV) 
URE 


(A, 
SUR 


0 


? 


NQAIN 
LV LOVE 
rr 


200 
fala l= 


CALCULATE VERTICAL AVERAGE PRESSURE 


C 


thd 
ac 
= 
Y) 
SY) 
Lil 
a 
Qa. 
-_ LL 
¥ ©) 
< <— 
ad aw 
mJ LJ 
CO > 
+ <{ 
ie 
<Z 
- Zz 
Soe Paks VY) 
- “NN => 
— = < 
= * —) 
> Ww = 
CQ it 
<{T—-_ 
+> Lh OL UL 
=D eDVDM& 


UOMO Kt Neate) 
PH CN tH aleHe ken 


= tezazsw 
OU DUS HYQRNIOE 
QON2xNGOAVOa 
mn o oo 
uy Oo ~mwo 
+ + +t 
oO 


EOS 


V(IeJSJ,K I-PBAR( I,J) 


OOORO0ONO 


Ooo 
Ooo 
MOPuy 


oT Xe TV yUM1, VM1 gUXG,UYGs VXGyVYG,UZG,VZG) 


o> 
~sJ 
> 
io 
Se 
Ni 
= 
co 
© 
a= 
-> 
ce 
i © 
=— < 
<f > 
me] Cc 
— J a2) 
a <I > 
CO —_ — 
a <— = 
<a © 
us mJ > 4 
ide _ ma 
ie ac * 
© > 
LL. ce es 
© — 
Y) ~ 
VY) Lis > 
t= ok he 
z <t os 
i al >< 
=< —_ = 
O W) ~ 
Qo. a FN ON 
=. q No Oa 
© Own o e 
WO —-Oxt << 
Ps OAZeEUO YO 
=> og UJ & ~ 
& Om Ond Dd 
= >OQOma «@ 
~ Ct —- — 
U. C—O OO 
Or WOaxrinia 
(Sj 2 eCDet) 
ZOD =“ S$-Oo>oO> 
O83 FRU 
mis “ase © 
E> OVUSIE—_ 
= com 
wed tl) Jia (ja 
Omg DWeOOW 
VY) fer > 
T 
wy 
OOOO 





oo 
Nout 


Newt 


(65 T=NS OSSETIA 


oO 
= 


- & = 


Next 
HU ome 


(6S T=n* CSF STON 


‘Oo 


SiH AD ae 
= 


paw DOw OO 


ee eOWre WwW AOOwW OW Qa 


((xt’9 


” noose Boe mont 
f— od (NU mt 
w—_ tho» om 


= 
YW wll led 


Ne 


-~ 
= 


000 © 
<rrr 
COW we 


Me 
6 0 
6 O 
6 OQ 
d N 


<I mt 4 4 
Onde 


LU OV OVO" 
- 


n 
(X 
(¢ 
(T 
17 
N) 


O 


(LWISSTWOSSASTWASN SE IW1S LWOS STdDSSASTVA SNS TV 


oO 


BSE 


L6E 


OT 


uy Om 
ey) 


uy 
‘0 


Qnuy 


OO 


104 


Cw fC IXY 


(LWOSSTWOSSASTWASNSTWN SE Td9S SLWOSSTVASA SIVA SN) 
S9 O1 OD (CT 


d3LS 3JWIL LXSN YOA SAIGVIYVA JONV4AD 
(MS9AANSOXNS ZOIMYSA 1179 
ALTIOIZA WIILYSA SILVINIIVI MYSA ANILNOYINS 


emint wee eNnmw W Dem © =O) ® =~ 


Ia>e> WOGOQN0D> ¢c«~ 
wQODT Me c© ONUOe YOu WO 


a a WZZZ REE F&F SMM Re OME OM 
POW) CBee OBZ FSF LOOaFt AQF AO 


wed weet Sti 


ra 


WwW 
amt 
Ov 


(9ZASIZNI DANS OK ASOANSOXNS 


TVA & X 
Sa 3 


Zu 
Oe > 
mel I 
ft J 
= a aad 
ad i) J 
OJ 
WY) 


TVNSALSX1§ 70% Vv9§ 
ALIQDOIZA WWINVZ 
LIDOWA IVIILYSA 


= 


TVAST VAS TWAS TWN SALSXLIENZASOZNIOANS DANS OXALOXN FAQVS 70% 120 
CTVASTVNFEAFENSALEXLIDZASIZNIIOANS DANS OXASOXNSAGV § 20 F140 
(TVASTVNSTWASTWNSALSXLIOZASOZNASOAANS SANS ah CN aS 
VV 

S 

S 


O Oo 
Oo Pf 
uy =u 


axis 
Seer 
J1LV1N 


NO 


I 
SLN3SNOdNOD SASN S 


® enw 


w~wXheJ J 
== IULOUOOOU NTO 


® eZOetHscHi=x 
Igoe ToOnovwoO 
a: 
2Z0f Se Je J 
BRU) aAugagod 


A 
3 
2 
Zz 
9 


uy 
=~ 
Oo 


OO 


Oo 





1.05 


JN3 
dOLS 
(9S THNSONSF FE TIIWIS) (OTS9)ASLIYUM 66€ 





INNILNOD 
(T-WH)Gx(OS-ISIF(CT-W SF ST) Ud=( SF ST) Ud 
COS OL U9 

(TWN) GadSt+( T=“ or ST) ad= C4 Sf ST) ad 
CS-—°2/CCNSF FETISH CT-N SFE 1) S)EdS 

TOT OL O9 (9° 29°H) AI 


wel 
QO O00OcCod0 
OOW NMNAaOAOUW 


YY) 


°Z/(T)Gxt(OS-( TSI 


0OO em 
a OQAtAle ee 
ut 


i-VWOII 


IW” 
{aw 
<tiual 
= UL ew 
OW e 
WY 


Y) 


Q3L93193N AYV SNOILVIYVA 3 
(6TIAS(9°TT’IE)S? 

; (6 jt 

(dd fa°s*v) 


w_—— DX 
j= ot awe 

2wMM wW 
O2Zz2Zz WY 
Yiwu << 
ots ww 
Saeme =O 
OMOA D 


zu 
(T 
T2 
$3 


s+ OOO 
Oo Oooo 
st UVASO 


Ooo 


J = 


106 





SLBROUTINE VERW(DD,UG;VG,W) 


iy 
Awe 
e «(f 
Od « 
vt ae) ae 
—~ & med 
Oso 
aN et 
weer ee eet CN] 
QWW~ 
OonN> z= 


aces. 
O000 
id ee 
NNW) 
222m 
UJ LU UL LL 
L=zz— 
a 
QAa00 


W SURFACE 


C 


=0 AND W BOTTOM=0 


s+ O- 
sO e& 
oe «& Il 
ede ij — 
Wm 
~x~Y eDD 
eMIZZ 
OoOo~7 emi 
OO] erdfef- 
SIN ZZ 
Oo0~w—WO 
QO FzBOO 


Ooo 
0 @ 


Ol 


oo 
Nest 


~ & 
NA 
i 


oo 


etn 


OO 
Oo 


0.0 


DO 300 K=1,8 


WIT 9J29) 


2*K 


K K 


vd bo oy 


—r ¢ 
wa eer CY 


it wi~— 
WO UL} = 


W(T gS eK D=WC 19S 2K-1)4+D0D (KK) *(CD+EF) 


TO 300 


GG 


DD(KK) ¥(CD+FEF ) 





IES AT 
PEN BOUNDARIES 


ELOCITIES A 
THE OPEN B 


-"O 


OO 


VGl2s1 Ki i72e 


CB=(VG(13,11,K)4+VG(1L49119K) 1/2. 
DE=(VG( 14,11 eK) tVG(15511,K))/2- 


2 
+ 
IF(KeEQel) GO TO 65 


GC TO 500 


LOG 


SLBROUTINE HVGRD(BB,CAs,DDyT Xe TY s UUs VV _UX gUY pVXeVV9UZs V2) 


ee caw =e 
DOD o& 
~ wn «SS 
© amt eal 
amie o 
= ee of 
AOOoOn 
ANIA a 
wwe et ee wea J 
Q> >>> 
aQ>TD> e 
oe 8 & Gem 
oe) oe eee, om 
= OO OV OW ony 
> a en ea 
—-<LOO0° 
OM) tl = tod 
~=a nn em & & 
O-OaAH6 
NUUNIAIAIOC 
ee SS ee ee SS gt 
<OD>x*<*<N 
O52 D>D 


Meee eee 
ON0000 
pang 4 04 9g Ped FG 
NIM NIN? 
2 Ae. 
UI LIS LL Li LLL 
ZLIAZAS 
ml | and (ed See) Fmd 
OAaAOaOO0o0 


N 
| 
Lt) 
= 
LN = 
Oe esQ> 
50 ¢ e~ 
Ol! AOO Il 
pin nha 
=> > IK J 
W>w~OOW 


5 
VL*¥RL)/EV 
DEFINE ZONAL GRADIENTS FOR 


INTERIOR 


C 


UXCI +1, J, KI=(UUCT 41,5 9K )-UUCI »J2K))/DX 
VX(I+1, 59K H(VVCI +l, Je KI-VV(II 959K) )/0X 


109 


LULU Ly 
Sa 
mee 


— —) =) 
he 
Zae 
OOO 


Ooo 


Ned 


FOR THE INTERIOR 


DEFINE MERIDIONAL GRADIENTS 


C 


>> 
ee) 
~N 
wx -_ 
ao & ) 
oe us 
oo = >< 
— <4 CQ 
— ow CQ ~ 
>> = << 
> <{ A-Oa 
1 | ~yeN~ 
om, i a om 
oe ” =——“ 
oo we <I Y~ ee 
= Liu QO - 
+ + “aN 
oe ws ews 
oo om <t aD: 
od == 
—— ” * > 
> _= > e> 
> Pd % NF 
oS ww we @uwr @ 
Qa it i mtv] | N 
Oo & mom CQ - aw jj — 
od eof we] 2 THeaq i) — jl 
nuh ww ce i mao 
ot YY Aa LULU YY or 
++ D> “| © 


ANM™MZZZWOO™” eM 
OOO eo eseianZ OO er o 
emt ward eer) tard band fee fe frm et SL et CN od 
OOO} S-ONOWOQOK KX 
QOOD>OVOCOO0O2D2>> 


VN 
OOo 
smi 


WO 


Kd) =-(2e*VV(20 559K) )/DX 


MLL 
——) 

“Za 
ound t=4 pang 
AF 
~Z2Z 
x<xOO 
mUOW 


OO 
Ooo 
LAS 


OUNDARIES 


— = — > 
=) > 
D+ > 
>erme 
HANNE N 

© e= Qu 
WNON 1 Od | 
= oe enews |] ~~ |] 
Occ ll — il — 
Cit Nomen Y 
Cer ey aL 
9 ee ert > 

OS8SadeiwtaZaa 
USO ew & @ arte 
ZO Pt tt tt ee b> 
— ee ee eee et 
LOO>>>>-OAY0 
WOoo7>D>> ON 


QO 
Ooo 
Of 
ho 


NTS AT THE NORTH AND SOUTH 


O 


é 


CALCULATE GRADIENTS AT THE OPEN B 


C 


(BB(K)-VV(15,10,K))/DY 
VV(1L,2,KI-CA(K)) /DY 


VY) 
— 
z 
> 
om | a, Ay 
Qa. Wey 
vu 
UL ww VY) 
SD COO al 
< Oa O 
cz ~™ veg 
Oo = - 
‘mee BD ey a= 
Y) vw OQ 
a & ae 
>O = Oo 
coc ~ o WO 
ome LL | 
On —~ > 
=O => w 
tL) J) => < 
m4 LL J = i CQ 
a> —_—a =< 
<I wh peed) = 
cwOWM | © 
COUILLI -_ pi fea} 
<_w ==) oo 
a) 4 ~— =o = 
fu”) > eo © 
ou - — = aaa 
ee 2 Pa >> © 
CW=AOoO >> @ 
Wt SHNOO ww a 
><_ eel ~ jl il Og 
aatetcn Al zx 


CLS Dm <TD med ND med SEI LL ULI << pee 
re Zmet ey | «© eDDD 
<< OO > Ye2BMZZ2zZuUIW 
AOmMOO “Ot @ etme 
A) ODN ed 4 0) CA) pt 24 bee free form <I tt 
C) <= mee ere ZILLI 
—+ COUNOZSNNQOOCA 
IMUOO>O~CDIS>OOVUID= 
Ore = Sm 
Ooo 
Main 


OVO WO 


LO 


—t + 
= @ 8 


ww NAIA 


OO << wor we 


Newdiee fl fl 
aexy—OO 
0D med oe oe LL 
— & & eDD 
oo 22 a ae ee 
OO & & eomtKiY 
LV P04 pg Pt few fe 
www ZZO 
OONNNQIOWaS 
OO>>DOVUCUW 


Ooo 
roy 





SLBROUTINE HORV (OT, PR »UXGyVXGeUYGe VYGyUZGeVZG eT Xe TY yUeVe UA VA) 


QAD22D2 Di 


Ze ZZ ZZ 
OQOOOOOU 
Od peng Pad md pees Dad Dad 
MMMM WY) 
2222222 
UL) OL LL LL 
> E> > i i a 
— me at el met 
Ooaoaaaaad 


MUN UY 


yet UY 


NONHO e ant 
©@ee OO €@ e 
OOOO li Oo 
0 0 de te rt 
her CK HS JAIN 
WW OO>mOW 


1 3 | | 
a 1) XN 
LEeoOCke 


Ae 


oo 
Dan 


om & 


| 
NZ“) = 


Oooo 
TON 


WOO 
O00 


KK=2*K 


PTY=—(L o/RNIXCPRIT +1, Jt ly KI-PROI ths Je KI +t PROT ep Utley KI-FPROIT eS eK) (26 


L*DY) 


Meee eee eee ee oer eee ee eee aia aaa ie aaa 


| > 


HOY 
1)/0 


oe eee Nau 


K ) 


(I,J 
*DT* UPTY+HD Y+VDY+WINY+UNL ) 
PRU tly Jy K)—PRU Ly Jp K)FPROUL +1, UtLy KI-PRUL, Jth eK) )/ (260 


N) *U 
)+2. 
)*VT 
N ) * ( 


cw aze 





(4° ST) YNKINY/TI) -=9LN 

LL/in=tH srs TI) yn 

( LAL+TNA+XNIM4+XQA +X GH4+X Ld Px LOx°Ze+( of SIIN=ALN 
(MoS TAXI NYS ZO) ETN 

(PF ALXIM-=XN 

(MHIGSOCTEN SP STIOZN-UHEK SE 1)9ZN) XC NYSAZ) ‘ 
xX 


( 
CNSR FE TIDAN-ONS THR S TIDANDFXO/CONSP SE TIOXN-( SKS T4412 9XN) D&C NY/HS ) 


it > 1 


A 
IM 
JA 
/ ( 
JH 


Oooo 
QoOoO 
NOS 


tele 





SUBROUTINE HDOGRD(0D,B0,CD,BByCAySySXySYySZ) 


Oo 
ee 


"ott 
xX 
aa 


OO 


CALCULATE ZONAL GRADIENT 


C 


JeK) 


a 


? 
9 


+1,53eK)-S(I1_.5,K))/0X 


~ 


1 


i 
x(3 


eOw <x) 
OatiN WW | 
ee em ! il 
aj ilo 
NL | ee 
a ioe | bee oe CETTE 
oe oe )DD 
ooe4H4) eZZz 
O90 + &AIni— 
PIN OO) ted = 
OOO << «KUO 
OOOMNMNMNOO 


oO Oo 


Oo Ooo 
o Nxt 


Lis 


CALCULATE MERIDIONAL GRADIENTS 


C 


eJt1l,K)-S(I1 959K) )/DY 


UV eg OL LU 
N + eNDD 
HQ00 Vee ZzZ 
QAHOIO © & wr 
Pm ed ad et et Le 
i ww 
YQOOO>>,+OO0 
NOOOWWMNOO 


2 Oo 
cr) Next 
= a= 


CALCULATE VERTICAL GRADIENTS 


C 


alot 
et NO 
=~ me & 
adie 
| 
el > 4 


ooowr 
eNO 3 
ANNAN 


UQO= 
OOoox 





(SCT sJeKI-SCI 9 Se K+L)) /DD(KKE1) 


+1) =(S(01 952K )-SC)/DOD(KK +1) 


© 

© 

iy 

© 

fom 

Oo © 

CO it © 
~ © 

= owl it 

0+6 ~ 


OO West LU 
GI eG «aD eDD 


LNONONGQO 
MMOMOMOUW 


Oo oOo 
Oc) Net 


MAI NIN 


CALCULATE MERIDIONAL GRADIENTS 
AT THE GRID POINTS IN THE VICINITY OF THE STRAITS 


C 
C 


OG 700 K=1,5 


CONTINUE 
U 


7100 


114 





SUBROUTINE ADVEC(DL»,BO,CD,UX, VY ,W2Z,BByCAySI ,ADV) 


THE BOUNDARY 


sKIV*E(STCI Ls KItSI(I+1,1,K) \l2e 


5°30 
1 


VY) 
~ 
-_ ras 
ro) = 
* © Oo 
~— oO. - 
ad 
& a 8 
od CO . 
AN =i rm ae 
N Lu = -— 
= = ( 
& Zz rug = 
—— oasg = ke 
OO & e = = 
oo LL A N + ©) 
mao le ~ ~ nad CO 
— wal = o=~ e ~~ => _dJ 
ov oe em 0 @ e e Ne ~e > LL 
et iO co AQIWZAA N N NN YN > > 
ww NA © ~ eX, ~ ~N ome, ~~ | 
- "YN> a aM se ety ey “MY > >> m4 
ey oo & © o | oa oo € ! YY J & - oe => r= O< 
one = a= “” ¥~w Qe Tt + omg LLY 
e000 - {-—— ee & oJ) et & (- - © 
ae © 2 >< omg ome One eno i= aE = Oz 
mMOLO 5 & —_) wate Owe QC) eS Tew ww =|) jC 
Ox > td > mde — Nae Se 2 LLJome 
ree eo & © ome /))U/) OW wer OOF OW) + =k 
N-—6 << <x+ + —+- mateo +t ++ i 
~ LALA my am ~— YY | $e e— 0. Xr& 
> www we tL OO +arYY © et e-+y Yn YOYY OD Oz 
OOax< Oo VCAN —I1 ee O- Qme ey er._ee i | Lc 
<ma> ota “aur bur ber P< < Oy 97) ered mest =t D> |! OO 
Zz NNA eee eQ><m™ @ wet en © O54! 1 Oem ® 
22z2zz OO WN MetArt OOH MOONY Wo. 
OD000 uy — wwe eo > | WK, © eet ewww 20. © OY) MLL em 
treme’ a ag LAG Ne 4 og pemeg Ora HE LL, fae pemeg {> 07) Oeeg PNY NY er eg OL. tt eee FE MD DD DO LI 
NAW) ston QIN OS wow OO Dee SE www ww YX OZZZ OM 
mee ee eee —+ OOOKX K Meee HO en mew wr > 1 I] Se I tet 
UL LULL LL oor 2D wNAesHNDDNNA LHR Y Me Nt oh tb fb Ww 
Lees won ut OO 1) wer nee were Sete er er eet ee ff OL re FO SZ 
ap Od eg XO + QOOx+tn nun nDdAwnuontnaranon ws >POUOUWME 
OQ00 oOaowW ZZ OOOYVOWLDPsVOULESA08>P> LN S IQ OO 
oO 
Oo VN emt 
or ooo 
=e, 
O QUO 


doi 





Kd+ST(LT-lsylsK) 7/26 
( 
S 


%* 1# YQ 
~x<Y~OMO e 

oD ex e& 
a—O ower © I 

one aaj eel i] i] 

ei i] oo] mY 

| me oe 

a ee | 

we OE OK ed ested LULU 
~<a Dette OD Det 
a> ~~ tl 


l 
? 
D 
ION OF THE ADVECTION AT THE SIDE WALLS 


Rox Moa oZZ— 

QIODIAHWOQ000Oa 

poe St pe Fe SIE KT OO 
Ooo 
oo 
NO 


OO 


ONS EXIST VELOCITY TANGENT 


TI 
ES 


OO 


Qi el 
=e Sw (/) S 
wet SJ) St er 
wr >) HE 
SN) ~~ | 
we | WIL 
o~<{ © 


e. 
Oma>elt- 0’ 
VOetelEem™! On 
ee i) « «a tim 
NAL « e=OOYrY - 
UUme CNA em Ym 
OD were wee “HLL red LL 
~~ >) > oleae &D 
OOr> er > ras raz 
Se0 nw leat 1 Naina 
LEN SPS SS wee LL LL ee em er er fee 
q{eoaseo moos 
O02 24aQ07 2QWOOO 
OO Re tee TOTO 


= © 


© © 
f uy 


116 


(uy) 


OO 


WIL 
OOx 


6 @ 
NG 
~N 
ay A 
et EE 
vam) oom] 
eel | 
NN 
~ “ & & 
ly Ay «cages A 
7 eee pant 
egjerf © & 
Lists 
Me ott 
CO ee oe 
| = a 
eo of/)Y) 
atl + + 
et ae? i, AED 
omg mt SS 
My) ea 
+ + 
Ay ES cool) aan) 
ww ee 
oe e fl 
= a 
a 
mnt (NI =t tog 
wow SU) 
dod Qed weer See” Lf"), 
WW) % O 
weet “ag omy, oe 
% H+ IO 
—~—l le 
ted Xi 
| «a «yD 
MS ett OD 
@ oe pal 
ae © oom 
e eFLnoO 
aAjt—t @ 
naw wet wae ws 
NNNNSI Wu 
=zzs e 
Wud as 
et NY eg eel me 
QQ dG) UL. 
IID Oe 


e@ ¢@ 
QIAN 
~ 
| le, A, 
EB) PB 
aml om) 
eee + 
NAY 
~~ NO & 
ae, iy go = 
om onl and 
ested o& & 
++ TK 
SM mt et 
eo @aw ww 
ad al ee ml 
oe H)Y) 
Han+ + 
3 Sew pt, SRD 
ee ae ee 
YN om 
++—1— 
ate et al ol 
<x oe 
oe fl 
aed eed =) 
eo ww 
ot Cle bt 
~~) 
el eee eee We” 
MAH Ft 
% HewX 
im iin Ge & 
SS eet ot 
© Gor— aml 
my om 
etl 
emt CC] eed ed OF) 
net eee eat woes 
NANIAIA 
x SK <x - -T ap) 
und ae 
at YOUN 
OvuMwoO 
AOaDAO 





ON 
NYUNLIY 


JNNILNID 

(MS TTSSTIAGV+ (44) IG/ ( ZEG-TEG)=(HS TTSSTIAIYV 
(NS TEST IAGV4+I HH) TOS (ZEN-TEND=HIHS TTS YT AIV 
(MS TSZVAQV4 (CH) 107 (290—-2dN )=(NS T £2 ADV 

(MS TST IAGVF (NH) TO/ (TIO-TdND=(4S TS TARY 


BC AUSSs | ei ST)IS)*O) weit’ 
“Z/AISIHINE TE CVIS DON'T 


eu) 


009 


Si 


Li 





SLBROUTINE DIFFO(DE,SXXsSYY,SZZ,_, ADD) 


OVO 


Ma 
it 
UL) Lis 
LL\ om 
LUVIN CG) U\ed 
NOOO eS e 
eeeme oD 
Tosi oo il 
“un von tw 
PLLA OU 
cCw<t<oOooaa 


ied 
OH 


one 
Ned = 
Wo 


ooo 


MetQm 


wood 
OOO 


weet ee (//) 
<>~— 
<> & 
nUY-= 
HH 
amy te VY 
<>~ 
OOwW 
™ SON 
om en, SJ 
MUL et 
#OORO 
Ot ed | 
a | | 
MK Oe NIN 
SH ONIN 


LS 


ACD(L ay JS eK) =XX4VYY+ZZ-ADD(I » Jy K) 


WLW 
> 
222 
(rat) (Re) f=) 
bem be 
Z2z 
VOW 
OOO 
coo 


SOO 
COAG ed 





= 

< 

~~ 
= we 
<I QV & 
~— =— 
S) = —_ 
wre CS 
ae Oo Oo 
f= =I <q w\ 
Mal Fe ON 
~w t= 
Oo ep 2 
OO t+ 
I< e 
oo N OO 
OOo om 
eo eo & Yam 
= lod oe o,) 
C2) wt pad “IV? 
= oo & a @ @ 
Fx =i am i 
WOO OF www (5 
| Za @ 
mz <I Ssteij—fe— 3 
Nie O OMNMWNOdO W 
AWA eee) Il ow tl 
LL Netcom 
Zae2Zx UN HM eOOO 
—SOOQOOWS et & wo ei 
fm png py =e CN] 3 w > ODD 


AW"NNHOOSOG e et e2ZZZ 
OZZS tQOOnm RO ete 
COS ELI LES LEY Se ed NI (0) er mee (/) fe wee fe fee fem 
VEE il aiqwv <dezZz 
2m OOO eH OF O0ORD 
NOOONMAOOANMHSONnVVUYO 


ooo°o 
Wooo 
NON et 


CONVECTIVE ADJUSMENT 


C 


<< 

<— 

x Ss 
Net Zu 
eee & 
wad od wed end 


tw 
= 


I ,J,K+tl ) 


0 
0 
0 
0 
K 
( 


1S 


B=STA(T 959K) 


+2 )*A)/(DD(KK)+DD(KK+2)) 


<~ = 

~ + 

ae ie 

oOo ~oo - 
WwW wor — 
~-_ -_ a 

O O00O = <t - 
= oe ae te) — 
+ =zOown < 

O agdaow ll | 
O F®OOU=- =~ 
— it = tt 


ZTOY~anN~yv~oroxod oo 

eOw 0 0 %® OO oO OOWO WWW 
fj—- = ODI DOVBOVDODINDDD>D 
Oo OWWe5wn ee ee ZF 222 
eO~— 6 CRO HOH ORK ORe 
CO bee HSS NC ee pee wes ee ee fe ff FF 
~~ Sereda <_« qt 2 222 
(LW iL Or Oe-OOODI0OUUY 
= OO SMNMNONIOHAYDIOQOOOUN 


© o Oo 98 8OoO 
wy oY O&O O99 
~~ OWT 





SLBROUTINE CVARCUN»,UNAs VN oVNAve SGNySGNAs UM, UNN 9 VMs VNN» SGMySGNN) 


Qov 

o oe 
oo 
=) = 

oo oe 
oo 
NIG 


ed 


<<a 
>> 


a 
iy 


~ & 


oo 
aad 

oo oe 
oo 
NAN 
ae 

=> >> eee 
e «DO 
a  & 
AQ ate 
e ented 
OOo ee 
oo ed 
e aWiC 
oo~~ 
NNIZ 
wwe 
—{ZzOO 
ZZM"wWM 
~~ ee 
oo om om 
am me 0 O 
OOD ew 
@ aed ad 
Oot 
aad © @ 
es ela 
OONN 
Nw~— 
~~Z2= 
Zed 
DDN 


zaza2ezaz 
O000 
at ee ed ed 
ONNMWY 
er 
LL LU tl La 
LL = 
ay) =O fame ot 
AOQOAaOa 


CHANGE HORIZANTAL VELOCITIES 


C 


oo 
Dan 
a f= & 
mimics 
ot tt 
~~" 


ooo 
et AIT) 


OOO 
COaa 


UM(I,JS,K) =UN( IT 9 J9K) 


UNN(IsJS 9K) =UNA(CT > Je K) 
VM(I gS e KI =HVNOIT 9 Jy K) 


/ 


VNN(I 99 J eK HVNACT p Se K) 


Ltt 
a= 
oa ae ee 
= 


[abe 


O00 
YOUO 


Oooo 


OooOe°o 
ON 


120 


CrhANGE SIGMA-TS 


C 


ade 
OO etN 


oe @ 


eae om 
Te 

ooo 

Fino 


GOO 
OOO 


SGM(I 9J 9K) =SGN(I 9 Jy K) 


SGNN(I,JsK)=SGNACI,J9K) 


LL LLL 
a> 
zzz 
= = 
Ee 


WUO 
OUOO 


Oooo 
OWNS 





Riot 2Or REPERENCES 


Bogdanova, A. K., 1961: The distribution of Mediter— 
ranean waters in the Black Sea, Okeanologiya, Wie, 
no. 6, p. 983-991; English trans., 1963, Deep Sea 


Research, v. 10, p. 665-672. 


Bogdanova, A. K., 1965: Seasonal fluctuations in the 
inflow and distribution of the Mediterranean waters of 
the Black Sea, Basic features of the geological struc— 
ture of the hydrologic regime and biology of the 
Mediterranean Sea, Edited by L. M. Fomin, Academy of 
Sciences, USSR, Moscow 1965. English trans., 1969 
Institute of Modern Languages, Washington, D.C. (AD 
686001). 


Bryan, K. and M. D. Cox, 1967: A numerical investiga-— 
DIGhCOtmene Oceanle general circulation. Tellus, 19, 
1, 54—80. 


Bryan, K. and M. D. Cox, 1968a: A non—linear model of 
an ocean driven by wind and differential heating: 


Part I. Journal of the Atmospheric Sciences, 25, 945- 
S67 ss 


Bryan, K. and M. D. Cox,,1968b: A non-linear model of 
an ocean driven by wind and differential heating: 

Part II. Journal of the Atmospheric Sciences, 25, 968— 
978. 


Defant, A., 1961: Physical Oceanography, Vol. l, 
Pergamon Press, 729 pp. 


Gunnerson, C. G. and E. Ozturgut, 1974: The Bosporus, 
In E. T. Degens and D. A. Ross, eds., The Black Sea — 


Geology, Chemistry and Biology. American Association 
of Petroleum Geologists, (Memoir 20), 1974. 

Haney, R. L., 1974: A numerical study of the response 
of an idealized ocean to large scale surface heat and 


momentum flux. Journal of Physical Oceanography, 4, 
145—167. — 


Levastu, T., S. Larson and K. Rabe, 1976: Numerical 
Models for Synoptic Analysis/Prediction of Currents and 
Temperature/Salinity Structure in the Oceans. Naval 
Research Reviews, June 1976. 


MZ he 





IEC. 


i. 


HZ 


Moller, L., 1928: Alfred Merz' Hydrographische Unter— 
suchungen im Bosporus und Dardanelenn. Veroffentlich— 


ungen Inst. Meereskunde an der Universitat Berlin, 
h@emholge wy Hett Fe Berlin 1928. 


Pektas, H., 1956: The influence of the Mediterranean 
Water on the Hydrography of the Black Sea: Technical 
Papers, 4th Meeting Fisheries Center, Istanbul, 1956. 


DEIyOtt, Pe, and O. Eigaz, 1946: The hydrography of 


the Bosporus an Introduction: Geo. Rev., v. 36, no. l, 
p. 44-60. 


le? 





PO. 


ii. 


INITIAL DISTRIBUTION LIST 
No. Copies 


Department of Oceanography, Code 68 .3 
Naval Postgraduate School 
Monterey, California 93940 


Prof. J. B. Wickham 3 
Department of Oceanography, Code 68 

Naval Postgraduate School 

Monterey, California 93940 


Prof. R. L. Haney ac 
Department of Meteorology, Code 63Hy 

Naval Postgraduate School 

Monterey, California 93940 


Deniz Kuvvetleri Komutanligi y 
Personel Egitim SB. Mudurlugu 
Ankara, Turkey 


Dz. Kuvvetleri Seyir ve Hidrografi Dairesi Bsk.3 
Cubuklu, Istanbul 


Turkey 

Istanbul Teknik Universitesl 1 
Taskisla, Istanbul 

Turkey 

Orta—Dogu Teknik Universitesi 1 


Ankara, Turkey 


Oceanographer of the Navy 1 
Hoffman Building No. 2 

200 Stovall Street 

Alexandria, Virginia 22332 


Office of Naval Research ih 
Code 480 
Arlington, Virginia 22217 


Dr. Robert E. Stevenson 1 
Scientific Liaison Office, ONR 

Scripps Institution of Oceanography 

La Jolla, California 92037 


Library, Code 3330 yz 


Naval Oceanographic Office 
Washington, D. C. 20373 


123 





iz 


i. 


mae 


nD. 


BEG 


ys 


18. 


neo. 


AUR 


2 


2s 


SIO Library 

University of California, San Diego 
Pao. BOx 2367 

La Jolla, California 92037 


Department of Oceanography Library 
University of Washington 
Seattle, Washington 98105 


Department of Oceanography Library 
Oregon State University 
Corvallis, Oregon 97331 


Commanding Officer 
Fleet Numerical Weather Central 
Monterey, California 93940 


Commanding Officer 

Navy Environmental Prediction 
Research Facility 

Monterey, California 93940 


Department of the Navy : 
Commander Oceanographic System Pacific 
Box 1390 | 

FPO San Francisco 96610 


Defense Documentation Center 
Cameron Station 
Alexandria, Virginia 22314 


Library (Code 0142) 
Naval Postgraduate School 
Monterey, California 93940 


Huseyin Yuce 

Zee t oo 

Deniz Harbokulu Egitim Uyesi 
Heybeliada, Istanbul 

Turkey 


Lcedr. J. F. Pfeiffer 

77 Adams Place 

Apt 407 

Quincy, Massachusetts 02169 


Lt. D. Henrickson 

Marine Sciences Branch 
Commandant (G—000) 

U. S. Coast Guard 
Washington, D. C. 20591 


124 











Thesis | | 166701 
| y825 Yuce . | 

mc. Numerical investigation 

of the dynamics of Sea of 
Marmara. 





thesY825 
umerical investigation of the dynamics 


DUDLEY KNOX LIBRARY 





